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We show that a recent claim that matter wave interferometers have a much higher sensitivity 
than laser interferometers for a comparable physical setup is unfounded. We point out where the 
mistake in the earlier analysis is made. We also disprove the claim that only a description based 
on the geodesic deviation equation can produce the correct physical result. The equations for the 
quantum dynamics of non-relativistic massive particles in a linearly perturbed spacetime derived 
here are useful for treating a wider class of related physical problems. A general discussion on the 
use of atom interferometers for the detection of gravitational waves is also provided. 

I. INTRODUCTION 

Chiao and Speliotopoulos proposed in a recent paper that gravitational wave detectors based on matter wave 
interferometry (MIGO) have sensitivity, under equal conditions, far superior to those based on laser interferometers 
such as LIGO. They argue that this is true especially for high-frequency gravitational waves and that it is better to 
use slow rather than fast atoms (see also Ref. 2] for related work). In a separate paper 3] Chiao and Speliotopoulos 
also argued that the usual quantum mechanical description for non-relativistic particles moving in the spacetime of 
a linearized gravitational wave is inconsistent. They suggest that the description in terms of geodesic 

coordinates (defined precisely in the next section) and that based on the geodesic deviation equation in the rigid frame 
are inequivalent. In this paper, we show that all those claims are incorrect. We also take this opportunity to address 
some general features of matter wave interferometry in comparison to light interferometry as applied to gravitational 
wave detection considerations. 

In the first part of this paper we show that the quantum mechanical description in terms of geodesic coordinates and 
that based on the geodesic deviation equation in the rigid frame are equivalent, contrary to Chiao and Speliotopoulos's 
claim. They just correspond to the use of two different coordinate systems and are related by a time-dependent unitary 
transformation. In the second part of this paper we compute the phase shift for an atom interferometer due to the 
passage of a gravitational wave in each of these two coordinate systems and show explicitly that the two results are the 
same (as one could have implicitly inferred from the equivalence of the two descriptions) . This contradicts the results 
of Chiao and Speliotopoulos in Ref. 0, where diflFerent expressions were obtained explicitly for the two approaches. 
We point out that their mistake originated from their procedure for computing the phase shift, resulting in different 
answers for the two coordinate systems. 

To help readers identify the source of confusion (perhaps at some risk of oversimplifying), we can explain their 
mistake as follows. Chiao and Speliotopoulos obtain the phase difference between the two ends of an interferometer 
arm by computing the total action for the classical trajectory with a given initial velocity that goes from one end 
of the arm to the other. In fact, they are interested in the change of the action due to the interaction with the 
gravitational wave, which is treated perturbatively. They keep the initial and final times fixed. However, when the 
perturbation is included, the position at the final time no longer coincides with the second end of the arm since the 
initial velocity and position also remain fixed. The correct result is obtained when the change of the initial time 
needed for preserving the boundary conditions is taken into account. 

When the correct analysis is employed, the great enhancement in the sensitivity of MIGO compared to LIGO found 
by Chiao and Speliotopoulos no longer exists, and it becomes clear that the sensitivity of MIGO is higher for faster 
atoms. Furthermore, despite the fact that the two cases correspond to slightly different physical situations (the atoms 
in MIGO are non-relativistic, whereas photons are highly relativistic), the expression for the phase shift in an atom 
interferometer that we obtain turns out to be formally equivalent to the general expression (for arbitrary gravitational 
wave frequencies) for a laser interferometer (if one replaces the laser wavelength by the de Broglie wavelength^ A and 
c with v) except for an additional term that accounts for the fact that the arms of MIGO are rigid, in contrast to 



Throughout this paper by de Broglie wavelength of an atom interferometer we mean the wavelength that free atoms with energy E 
would have in the absence of gravitational waves, where E is the constant energy at which the atoms are emitted by the interferometer 
source. 
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laser interferometers such as LIGO, whose mirrors are freely suspended. (The implications of having a MIGO with 
freely suspended mirrors are discussed in Sec. IVjl. 

Whereas LIGO operates in the low frequency regime (transit time smaller than the gravitational wave period) , for 
most relevant situations MIGO would operate in the high frequency regime. In that regime, the term in the expression 
for the phase shift due to having rigid arms dominates. That term is given by i/A, where L is the length of the arm, 
times an oscillatory factor of order one which is a harmonic function of the time times lo for a gravitational wave with 
angular frequency to. Therefore, under equal conditions, i.e., the same values for the wavelength, the length of the 
interferometer arms and the flux, the sensitivity of MIGO would be comparable to that of LIGO. 

Going beyond the correction of mistakes in previous analysis on the same problem, one can still ask whether there 
is a feasible range of parameters (wavelength, length of the interferometer arms, and flux) for an atom interferometer 
such that its sensitivity for the detection of gravitational waves is better than that of LIGO. We will briefly discuss 
this point in Sec. With these discussions as a start, it is our hope that some interest in both the gravitational 
physics and the AMO (atomic, molecular and optical physics) communities will be generated to explore in broader 
terms the scientific potential of atom interferometry toward precision measurements related to gravity effects. 

The paper is organized as follows. In Sec.|nlwe derive the geodesic equations describing the motion of free-falling 
particles in two different coordinate systems (the geodesic versus the rigid), and show how they are related to each 
other by a time-dependent local transformation. In Sec. Illll we write down the action for the non-relativistic motion 
of a massive particle in the spacctime of a linearized gravitational wave and derive the Hamiltonian in terms of the 
canonical variables for these two coordinate systems. In the quantum theory obtained from canonical quantization, 
we write down the Schrodinger equations for the two coordinate systems used and show that they produce unitarily 
related dynamics. This disproves the claim that only the rigid coordinates give the correct physical results. In Sec. II VI 
we derive the phase difference between the ends of an interferometer arm including the effect due to gravitational 
waves from an analysis of the Hamilton-Jacobi equation, which corresponds to the lowest order approximation for 
the Schrodinger equation in a JWKB semiclassical expansion, with the action expanded in powers of the metric 
perturbation. We show explicitly that the results for the phase differences calculated in these two coordinate systems 
coincide and identify where the error was made in earlier works. In Sec. we discuss the physical implications of the 
result obtained in Sec. II VI as well as some general aspects about the use of atom interferometers for the detection of 
gravitational waves. Some technical details about an alternative way to compute the phase differences are included in 
Appendix 1X1 They are useful in order to compare with Chiao and Speliotopoulos's approach. Finally, in Appendix FBI 
we explain why, contrary to the claim made in Ref. 10], the effects due to the reflection off the mirrors cancel out in 
the final result for the phase shift. 

Throughout this paper relativistic units with c — 1 are used. Indices are raised and lowered using the background 
Minkowski metric ry^j^. Greek indices correspond to spacetime indices running from to 3, whereas Latin indices 
correspond to spatial indices running from 1 to 3. The Einstein summation convention over repeated indices is 
employed. 



II. GEODESIC VERSUS RIGID COORDINATES FOR A LINEARIZED GRAVITATIONAL WAVE 

In the transverse-traceless gauge, the metric for a linearized gravitational wave propagating along the z direction 
is given by 

ds^ = -dt^ + {5ij + hij{t - z))dX^dX\ (1) 

where correspond to three spatial coordinates (in particular z = X^) and hij is a tensor-valued function with 
arbitrary dependence on < — z that satisfies the conditions /i^s = and h\=Q |^ . 

The geodesic equations describing the motion of free-falling particles in this coordinate system are 

d'^X' -^dtdX^ -^dzdX^ 1^. dX^dX\^.,^^ 

d^t 1; dX'dX^ ^,,2, 

where r is the proper time and the overdot denotes a derivative with respect to the inertial time t. The initial 
conditions must satisfy the normalization condition {dt/drY — i^ij + hij){dX'^ / dT){dX^ / dr) = 1, which is then 
preserved up to linear order in hij by Eqs. l|3l-lI2l- For a non-relativistic particle we have (dt/dr) = 1 + O(w^) with 
= {dXydt){dX,/dt), and Eq. 101 becomes 

d'^X' ■ dXi 



h)— + 0{hi^yh.,). (4) 



dt^ ~ J dt ' ^^"y'" 



3 



Since worldlines with X^{t) = constant are solutions of the non-relativistic equation (0J, or even Eq. we will use 
the term geodesic coordinates for this coordinate system (commonly referred as transverse-traceless coordinates). 
Alternatively, one can work with a new coordinate system related to (i,^*) by a transformation 



t = t, (5) 
1. 
2' 



^X^ + -h)X\ (6) 



which is generated by the vector field £^^{t,X^) ~ — (l/2)?7^'/iij(t — z)X^ . The metric perturbation in this new 
coordinate system becomes /i^^ = + d^S^^ + d^S^f^ and hence the full metric is given by 

ds^ ^ -dt^ - hij{t - z)x'-dtdx^ + Si^dx^dx^ + hij{t - z)x'dzdx^ + 0{h'^j), (7) 

where terms of quadratic or higher order in hij have been neglected. It should be noted that this coordinate system 
is only consistent with the linearized treatment of the metric perturbation as long as hij{t)x^ is of the same order 
as hij. Therefore, if we consider a gravitational wave of wavelength AgWj we should restrict our consideration to a 
region such that < Agw for any i. 

The geodesic equations in the {t,x^) coordinates are 

d'^x' I--, , f dtV ,dtdz 1..^ , f dzV 1 ^ . dx'dxi „ ^ 



dT^ 2 ' \dT) 2^ drdT 2^ \dTj 2 '' ''dTdT 

d^t 1; dx'^ dx^ ^.,9n 

Proceeding analogously to the previous case, one obtains the following equation for a non-relativistic particle: 

d^x' Iv,- . Iv, .dz 



From Eqs. 11U|) or Q one can see that in the non-relativistic regime it is possible to have solutions of the geodesic 
equations in both coordinate systems with z{t) = constant. For simplicity, from now on we will restrict ourselves to 
the case in which the motion of the non-relativistic particle is confined to the z = plane. Eq. (|10|l reduces then to 

^ = lhy+0{hl,v'h,,). (11) 

Since one can verify that the physical distance between worldlines with (t) = constant remains constant in time up 
to quadratic order in hij, we will use the term rigid coordinates for this coordinate system. 

Eq. (|ll|l can also be derived from the geodesic deviation equation which governs the separation of nearby 
geodesies and reduces for rigid coordinates to 



d^x' 



i?'oo,^^ (12) 



dt^ -^^00, 

by taking into account that the Riemann tensor to linear order in the metric perturbations is given simply by 

1 
2 



i?'oo, = :,h) . (13) 



However, this is not a peculiar feature due to using rigid coordinates: Eq. Q, for instance, can also be derived from 
the geodesic deviation equation expressed in terms of geodesic coordinates. 



III. DESCRIPTION OF A NON-RELATIVISTIC PARTICLE PROPAGATING ON A LINEARIZED 

GRAVITATIONAL WAVE 

The action, from which the geodesic equations can be derived, for a particle of mass m propagating on a spacetime 
characterized by a metric gp,i, is given by 



S[y'{j)] - -rn I dr^-m I djJ ~g^p{yf^{j))^^, (14) 
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where 7 is an arbitrary parameter for the particle's worldhne. Here we used the notation y'^ for the spacetime 
coordinates in a general coordinate system to distinguish them from the particular coordinate systems used to describe 
a linearized gravitational wave introduced in the previous section. The action for a non-relativistic particle propagating 
on a linearly perturbed Minkowski spacetime with metric g^v{x) = Tjf^u+h^y(%j), can be derived by taking the particular 
parametrization 7 = j/*^ = i and neglecting terms of order hfj^ or higher. In particular, if /iqo = (as will be the 
case in the situations addressed in this paper) the non-relativistic action becomes 

Sly'm . / ,tL (m.m - /■« I (f f + 2/;..*, „<)f + i.,('.v')'if - 2) ^ (15) 

The last term in the integrand can be rewritten as a time derivative oimt and does not affect the equations of motion 
derived from the action. 

We now consider the spacetime of a linearized gravitational wave in the two coordinate systems introduced in the 
previous section. Recall that we restrict our attention to the case in which the motion is confined to the z = plane. 

A. Geodesic coordinates 

When geodesic coordinates are used, the metric is given by Eq. |^ and the Lagrangian in Eq. ((Tsl) reduces to 

L{X\X') - ^ (X'X, + h,j{t)X'X^ - 2) . (16) 

The equation of motion derived from this Lagrangian coincides, as it should, with Eq. which was obtained directly 
from the geodesic equation. 

One can easily make the transition to the canonical formalism. With the conjugate momenta Pi = rn(Xi + hij{t)X^) 
one obtains the following result for the Hamiltonian up to linear order in the metric perturbation: 

H^{X\P,) = ^P,P' - ^h,j(t)P'P^ +m. (17) 
2m 2m 

B. Rigid coordinates 

Alternatively, if the rigid coordinates are employed, the metric is given by Eq. (|7|, and the Lagrangian for the 
motion of a non-relativistic particle becomes 

L{x\x^) = ^ [x'x, - h,j{t)x^x' - 2) . (18) 

The equation of motion derived from this Lagrangian coincides with that obtained from the geodesic equation and 
given by Eq. l(TT|) . 

Again, one can easily make the transition to the canonical formalism. The conjugate momenta are pi = m{xi — 
(l/2)hij{t)x^) and the Hamiltonian up to linear order in the metric perturbation is 

H,{x\pi) - -^p.f + h,j{t)x'p> + m. (19) 
2m 2 

It should be noted that the last term in Eqs. H17|) and (|19|) simply corresponds to the rest energy of the massive 
particle and can be absorbed by redefining the origin of energies. When doing so, at the quantum level one should 
also redefine the state vector |V'(i)) to eKp(imt)\'ip{t)) . From now on we will presume that such redefinitions have been 
made. 



C. Equivalence between the two descriptions 

The two Lagrangians are equivalent, up to linear order in the metric perturbation, under the time dependent 
coordinate transformation 



X' ^ x' = X' + -h){t)X^ . 

2 



(20) 
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This change of coordinates induces the following transformation for the momenta up to linear order in the metric 
perturbation: 

P^^p^ = P^^lh,,{t)P\ (21) 

which together with Eq. (|20|l constitutes a time dependent canonical transformation, i.e., it leaves the symplectic 
structure of phase space (or, equivalently, the Poisson brackets for the position and momentum variables) invariant 
at every instant of time, up to that order. Under this canonical transformation (X* — > a;*, — > Pi) the Hamiltonian 
transforms, up to linear order in hij, as 

Hg{X\P,) ^ H,{x\p,) = H^{X\Pi) + h,,{t)X'P^, (22) 

where the subscripts g and r denote quantities in geodesic and rigid coordinates respectively, but PiX'^{X^ , Pj) — 
H{X^,Pi) remains invariant. 

When quantizing the theory by promoting the position and momentum variables to operators satisfying the canonical 
commutation relations, the previous canonical transformation corresponds to an analogous transformation for the 
position and momentum operators that preserves their commutation relations. According to the Stone- Von Neumann 
theorem (see, for instance, Ref. 12] and references therein), this is sufficient to guarantee that the two representations 
are unitarily equivalent, i.e., connected by a unitary transformation, which will be time dependent in this case. 

In fact, the unitary equivalence of the descriptions can be shown explicitly. Let us start with the Schrodinger 
equation in terms of the geodesic coordinates 

where we have introduced the position basis^ l2/*)gj characterized by Xi\y'-)g ~ y-'W)g and g{y^\y"')g — S'^^^y^ — y"-). 
The action of the momentum operator in this basis is given by g{y^\P-'\y'^)g — —ih{d / dy^)6'^'^\y^ — y'^)- Similarly, one 
can consider the position basis in rigid coordinates |2/')r, which is characterized by x^\y^)-c = y-'\y^)x and r{y^\y'^)r ~ 
^(3)(J/^ _ yHy Taking into account that = X' + ^h){t)X^ , one finds that the position basis in the two descriptions 
are related up to linear order in hij as follows: 

\f),^\f-{l/2)mt)y^),^U{t)\f)„ (24) 

where the linear operator U{t) relating the two bases is unitary up to linear order in hij since it preserves both 
orthogonality and normalization up to that order, as can be seen from Eq. H24|l together with the fact that the value 
of the Jacobian for the coordinate transformation in Eq. (|20|l is one plus terms quadratic in hij , which is a consequence 
of hij being traceless. The same unitary operator also relates the position and momentum operators within the two 
descriptions in the following way: a;* = U{t)X^U^^{t) and pi = U{t)PiU~^{t). One can check that it is given by 
U{t) = 1 + {i/2)hij{t)X^P^ + 0{h1j). It should be noted that the ordering of the operator product X'^Pj appearing 

in both U{t) and is irrelevant because it is always contracted with hij and hij{t)[X^ , P^] vanishes because hij is 
traceless. 

Furthermore, the Schrodinger equation in terms of rigid coordinates, 

ihj^Mit)) ^ .{y'\Ms^\Pi)m))- (25) 

is equivalent, under the time-dependent unitary transformation generated by U(t), to the Schrodinger equation in 
terms of geodesic coordinates, given by Eq. H23|) . This can be shown as follows. First, we rewrite r(y'|'0(O) as 
/ d^y'r(2/'|y'*)gg(2/"|'/'(i)) on the left-hand side of Eq. 123), which becomes 

I '^V(^r(y%'^>g) g{y'^m)+^f^ I d'y'Ayy')gj^g{y''\m)- (26) 



^ A mathematically more rigorous treatment sliould employ strictly normalizable states, but identical conclusions would be obtained if a 
basis of normalizable wave-packets were used. In addition, we should keep in mind that, when using rigid coordinates, the perturbative 
expansion in hij is only valid for regions smaller than the wavelength of the gravitational wave, as mentioned in Sec. in] Therefore, 
throughout this subsection it is implicitly assumed that the Hilbert space is restricted to states localized within a region smaller than 
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Next, we use the operator version of Eq. (|22() to substitute Hi-{x'^,pi) into the right-hand side of Eq. (|25|l . which 
becomes 

J d^y\{yW^^{y'^H^{X\P.,)m)) + ^h^kit) J d^y\W%s{y''\X^ P'lm)- (27) 

The last term in this expression is exactly canceled by the first term in expression H26() since ih{d/dt) ,-{y''\y''^)g — 
{ih/2)hjk{t){d/dy')5{y - {l/2)h{t)y - y') = {l/2)hjk{t)g{y''\Xi P^\i^{t)) , which follows from Eq. The whole 

equation becomes then 

^n J d^y',{f\yX-^^siy"\m) = / d^y\W), ,{y'^\H,{X\ ■ (28) 

Finally, multiplying both members by g(y"'|2/*)r and integrating over one recovers Eq. (|23|l . Therefore, the time- 
dependent unitary transformation generated by U{t) accounts for the different time evolution generated by the cor- 
responding Schrodinger equations for the two different Hamiltonians H^{X'^,Pi) and Hi.(x'^,pi). 

To summarize, it is clear that, contrary to the claims made in Refs. 3, 13], the quantum mechanical description in 
terms of geodesic coordinates is equivalent to that based on the geodesic deviation equation, which simply corresponds 
to using rigid coordinates, and there is no inconsistency with the usual approach, which employs geodesic coordinates. 



IV. PHASE SHIFT IN ATOM INTERFEROMETERS CAUSED BY GRAVITATIONAL WAVES 

In this section we will analyze the horizontal MIGO configuration proposed in Ref. 1], in which the Earth's 
gravitational filed plays no major role and its effect is neglected. It is slightly simpler than the vertical configuration 
also proposed there, but it already contains all the issues that need to be addressed. Throughout this section we will 
consider interferometers with rigid arms, as proposed in Ref. 0. Additional discussion on the implications of having 
an atom interferometer with freely suspended mirrors is provided in Sec. 



A. General expression for non-relativistic atoms 



The phase shift due to the passage of a linearized gravitational wave through an atom interferometer can be 
computed using the geometric optics approximation (the lowest order JWKB approximation) as long as the wavelength 
of the gravitational wave and the arms of the interferometer are much longer than the de Broglie wavelength of the 
atoms. Under this approximation one only needs to solve the Hamilton-Jacobi equation 



dS_ ^ f dS 

dt \ dy^ 



^^—+H[—.y\h,,{t)], (29) 



perturbatively with respect to the metric perturbation. Here we have used a generic notation y*, which can correspond 
to either X^ or x*. Expanding the action in powers of the metric perturbation as S{y^ , t) — Sa{y^, t) + Si{y'^,t; hij] + 
0{h^j), we obtain the following set of equations from Eq. H29I) : 



55*0 


1 


fdSo' 


dt 


^ 2^ 


\dy\ 


dSi 






'df " 


m \ 


.dy' J 



dSo 



6'^ (30) 

(31) 



where H^^"^ {dSo/dy^ ,y^ , hij(t)) corresponds to the interaction term of the Hamiltonian, which is linear in the metric 
perturbation. 

Eq. H30|) can be solved by the variable separation method and the result, which corresponds to the phase of an 
energy eigenfunction in the semiclassical approximation, is 

Soiy\ t) = -Et + 2TTjn,f + h^bo, (32) 

where E is an integration constant that corresponds to the energy of the particle, A = h{2mE)~^/'^ is the corresponding 
de Broglie wavelength, is a unit vector in the direction of wave propagation and <f>Q is the initial phase at t = 
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for any point on the plane defined by riiy^ = 0. The expression for So{y'^,t) is the same when using either geodesic 
or rigid coordinates since the part of the Hamiltonian independent of the metric perturbation has the same form in 
both descriptions. However, one has to take into account that in general the trajectories of the different elements of 
the interferometer, such as mirrors and beam-splitters, will be different in each coordinate system. 
Substituting Eq. into Eq. we get 

^ = -H('^(^jn,,y\h,,{t)y (33) 

where dSi/dt = dSi/dt + v^dSi/dy^ with = {h/mX)n^. Eq. 1)33(1 has the same form for both the description using 
geodesic coordinates and that employing rigid coordinates; however, the explicit expressions for If^^^ differ. 

The phase difference, up to linear order in the metric perturbation, between the two ends of the arm of an atomic 
interferometer of length L for non-relativistic atoms of energy E propagating along the arm is A(j){t) = ASq {t)/h + 
ASi{t)/h with 

ASo{t) = So{yh,t) - So{y\,t- L/v) = ^n, {y],{t) - y\(t - Ljv)) + (34) 

A V 

where y\ and are the positions of the two ends of the arms, which will be different for each coordinate system, 
and 

i 

ASx{t) = Sx{y'j,,t)-Sx{y\,t^Llv)^- J dt' H^^^ (jn,,y\^^{t'),h,,(t')^ , (35) 

t-L/v 

where y\Q-^{t') — y\ + {vt')n^ and coincides with the classical trajectory for a particle of mass m and energy E in the 

absence of the gravitational wave. Furthermore, we can take the value for y\ in either coordinate system because the 
difference corresponds to terms that are at least linear in the metric perturbation and would yield a contribution to 
Eq. H35|l of quadratic or higher order in the metric perturbation. 

Eqs. l|M)) and agree with the general results obtained in Refs. 0,0. 



B. Explicit result using geodesic coordinates 



If we assume that the arm of the interferometer is perfectly rigid, i.e., that the physical distance between the two 
ends remains constant in time, their positions in rigid coordinates, x\ and Xg, are constant'^. Their positions in 
geodesic coordinates can be easily obtained, up to linear order in hij, by inverting Eq. 



X\{t - L/v) = x\ - -h){t - L/v) x\ 

XW) = x^B-lhmx's. 



Thus, in geodesic coordinates, ASo{t) is given by 

ASoit) 



h 








hL 













Xa) - -n 



-K,it')yl^^{t') 



\h.,j{t')yl„^{t') 
t'=t 

+ 

t'=t~L/v 



t'=t 

t'=t-L/v 



(36) 



L 
E- 

V 



(37) 



where y^/^Q-^it) = and y-j^^-^it^ L/v) = x\, and wc took into account that ni{x'^g ~x^a) = ^- Substituting H''^\ which 
is given in this case by the second term on the right-hand side of Eq. H17|l , into Eq. H35|l , one obtains the following 



^ This is true provided that we take the origin of coordinates at the center of mass of the interferometer and it follows a geodesic. In fact, 
if the gravitational wave propagates perpendicular to the interferometer, it is sufficient that its center of mass can move freely along 
any transverse direction, i.e., on the plane determined by the two arms. 
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result for ASi{t): 



(38) 



t-L/v 

Therefore, the total phase difference between the two ends of the interferometer arm is 

/ 



A<^(t) = 2tt 



L 

A " T 



V 



n t'=t 
t'=t-L/v 



V 

2A 



dt'h,j{t')n'n^ 



E L 

h V 



(39) 



t-L/v 



C. Explicit result using rigid coordinates 



by 



When rigid coordinates are employed the positions of the arm ends are fixed at and xb, whence AS'o(i) is given 



ASo{t) = -n, (x'b 



L hL L 
E- = — + E-. 



(40) 



Substituting H^^\ which corresponds to the second term on the right-hand side of Eq. (|19|l . into Eq. H35|l . the result 
for ASi{t) is 



ASi{t) 



h 
2A 



t-L/v 



hv 
2A 



t-L/v 



2A 



t'=t-L/v 



(41) 



where we have integrated by parts and taken into account that y^Qj(f') = vn^ . Finally, the total phase difference is 

+ ^ 'j dt'K,{f)n^nn +1^, (42) 



t'=t 

t'=t-L/v 



t-L/v 



which coincides with Eq. 

The first term on the right-hand side of Eq. H42|) . which is independent of hij and does not change in time, 
simply corresponds to the phase shift in the absence of gravitational waves, whose effect is contained entirely in the 
remaining two terms. If we assume, for simplicity, that the first end of the arm coincides with the center of mass of 
the interferometer with rigid arms, and the interferometer is suspended in such a way that its center of mass can move 
freely along the transverse directions (i.e., along any direction on the plane determined by the two arms), we can 
take XA = 0. The second term on the right-hand side of Eq. (|42|1 is then proportional to {L/\)n^n^hij{t). Finally, if 
we consider the particular case of a monochromatic gravitational wave of angular frequency lo, the last term becomes 
proportional to v/Xlo times a dimensionless oscillatory factor, which is a product of a harmonic function of cut and a 
harmonic function oiujL/v. 



D. Phase shift 



In the previous subsections we obtained the phase difference between the two ends of an interferometer arm. Here 
we explain how to use those results to compute the phase shift between the two interfering waves in an interferometer. 

First, one divides the trajectory followed by each interfering wave in stages that describe free propagation between 
the different elements of the interferometer: source, beam-splitters, mirrors and detector. Next, one takes one of the 
interfering waves and adds the phase differences for every stage (using Eq. (|39|l 'l to obtain the difference between 
the phase at the detector location and at a given observation time tf, and the phase at the source location and time 
ti — Lj/v, where Lj is the total length (in the absence of gravitational waves) of the path followed by that 
wave. Similarly, one can get the analogous phase difference for the second interfering wave between the observation 
time t{ and if — J2j -^j/"' where J2j is the total length of the path followed by that wave. Finally, one needs the 
phase difference, for a fixed position at the source location, between times tf — J2j Lj/v and if — J2j -^j/^- adding 
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the latter to the phase difference for the second wave and comparing with the phase difference for the first wave, one 
obtains the phase shift between the two interfering waves. 

We now turn to the computation of the phase difference between times t[ — '^jLj/v and if — '^jL'^/v at the 
source location. For simplicity we will consider the situation in which the source is located at the center of mass 
of the interferometer and take into account the remarks made in footnote 2. It is then natural to assume that the 
atoms are emitted at constant velocity in the rigid frame if the size of the source is much smaller than the arms of 
the interferometer. This is because the term that accounts for the interaction with the gravitational wave in rigid 
coordinates, the second term on the right-hand side of Eq. Ijiyf) . is proportional to the distance from the center of mass 
and can be neglected provided that the size of the source is small enough. Furthermore, under those conditions it is 
also a good approximation to assume that atoms emitted at constant velocity in the rigid frame correspond to atoms 
emitted at constant total energy E since the interaction term is negligible. Hence, from Eq. H29|l one can conclude 
that the phase difference between times tf — Lj and if — L'^ /u at the source location is given by 




(43) 



When we add this to the phase difference for the second wave and subtract the result from the phase difference for 
the first wave to obtain the phase shift between the two interfering waves, the contribution from the last term on 
the right-hand side of Eqs. I|39|l or H42|l . when added for all the stages of the path followed by each wave, cancels the 
right-hand side of Eq. 14311 exactly. Thus, when computing the phase shift, it is enough to add the phase differences 
for each stage keeping only the first three terms on right-hand side of Eqs. (|39|l or l|42|l . Note that the whole argument 
also applies when using geodesic coordinates, even though the atoms are not emitted at constant velocity in those 
coordinates, because the difference between the Hamiltonians in rigid coordinates and geodesic coordinates is of the 
same form as the interaction term in rigid coordinates (see Eq. I|22|l ). Therefore, whenever that interaction term can 
be neglected for the reasons explained above, one can assume that the atoms are emitted at constant energy in terms 
of the geodesic coordinates as well. 

To illustrate how the phase shift is computed with a particular example, let us consider an interferometer with a 
Michelson-type configuration with arms of length L and L', and a negligible distance from the atom source to the 
beam-splitter as compared to the arm lengths. One needs to add the two contributions that correspond to the two 
trips in opposite directions along one of the arms using just the first three terms on the right-hand side of Eqs. H39(l 
or (|42|l . proceed analogously for the second arm, and finally subtract both results. Assuming that the source and 
beam-splitter are located at the center of mass and that the interferometer is freely suspended as specified in footnote 
2, the result for the phase shift is 
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(44) 



where fc' and P are the unit vectors parallel to each arm. For the particular case in which the two arms are perpen- 



dicular and have equal length and hij{t) is a harmonic function, i.e. 

for the phase shift becomes 



hij{t) = hij sin^Lut + ip), the previous expression 
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(45) 



E. Discrepancy with previous results 

There is a discrepancy between the result for the phase difference obtained in Sees. IIV Bl and IIV CI and that 
obtained in Ref. 0, which contains an extra term proportional to V^uj/Xv times a dimensionless oscillatory factor 
that dominates in the regime considered there. In this subsection we explain the relationship between our approach 
and that of Refs. for computing the phase difference between the ends of an interferometer arm. We also point 

out a mistake in their analysis which explains the discrepancy with our result. 
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The basic quantities in our computation and in that of Ref. [lj,|9( are shghtly different. Here we have considered the 
solution of the Hamilton- Jacobi equation at fixed energy to lowest order, whereas in Ref. the action for a classical 
trajectory that goes from xa to xb in a given time was employed. Although these two objects share some common 
properties (they both satisfy the Hamilton- Jacobi equation H29I) ) they have different meanings. This point can be 
easily illustrated by considering the example of a onc-dimcnsional time-independent Hamiltonian. The analog of the 
object that we considered is a solution of the Hamilton- Jacobi equation of fixed energy E, i.e., a function Ssixji) 
which satisfies the Hamilton- Jacobi equation and such that dSE{x,t)/dt = —E. On the other hand, the analog of 
the object considered in Ref. Xn M the function S{x,x',t) which results from evaluating the classical action for a 
classical trajectory that goes from x' to a; in a period of time t. It can also be obtained from Se{x, t) by choosing the 
energy of the classical trajectory characterized by x' , x and t: S{x,x' ,t) = Se{x,x' ,t){x,t)- Although one can show 
that S(x,x',t) also satisfies the Hamilton- Jacobi equation, SEix,t) and S{x,x',t) are not equivalent as functions of x 
and t. Furthermore, whereas SE{x,t) can be identified with the exponent of the lowest order JWKB approximation 
for an eigenfunction of energy E: 

V'£;(x,t) ~e*•^^(^■*)/^ (46) 
S{x,x',t) corresponds to the exponent of the lowest order semiclassical approximation of the propagator: 

K{x,x',t) = 5]^i^(x,t)^|;(x',0) ~ e^^^(-'-''*)/^ (47) 

E 

where the sum is over all the energy eigenvalues, including possible degeneracies, associated with the Hamiltonian 
that one is considering. Both iIie{x, t) and K{x, x' , t) satisfy the same time dependent Schrodinger equation, but they 
are again different functions of x and t. Since the physical situation that we are interested in is the interference of 
atoms emitted with a fixed energy (up to a good approximation) propagating along the arms of the interferometer, 
TpE{x,t) (and SE{x,t) in the semiclassical approximation) is the relevant quantity that should be employed to obtain 
the phase difference directly. 

Nevertheless, one can still make use of S{x,x',t) to compute phase differences for an interferometer. This follows 
from the observation that by evaluating the action for a classical trajectory of energy E that goes from x to x' one 
gets S{x,x' ,T) — Se{x' ,i) — SE{x,t — T), where T is the time needed to go from x to x' . The phase SE{xf,t) at 
the detector location for one of the interfering waves in an interferometer computed in this way is given by the action 
for the whole trajectory that corresponds to the path followed by that wave in the interferometer plus SE{xi,t — T), 
where Xi is the position of the atom source and T is the time the whole trajectory spends in the interferometer. The 
phase shift between the two interfering waves is then given by the difference of the actions for the classical trajectories 
that correspond to the two paths followed by the two interfering waves plus SE{xi,t — T2) — SE{xi, t — T{), where Ti 
and T2 are the times spent by the classical trajectories in the interferometer. For a time-independent Hamiltonian 
the last contribution simply corresponds to E{T2 — Ti). 

A similar procedure can be employed when a time-dependent perturbation is included in the Hamiltonian. In 
this case one needs to evaluate perturbatively the action for the classical trajectory associated with each interfering 
wave, which is also determined perturbatively. When computing the phase difference between the two ends of an 
interferometer arm, one should consider solutions of the perturbed equations of motion that go from one end of the 
arm to the other. Furthermore, it should also be taken into account that the kinematic momentum is reflected on the 
mirrors and, when using rigid coordinates, the atoms are emitted with constant velocity, as explained in Sec. II V PI 
For a given observation time at the detector, the previous conditions completely determine the trajectory. That means 
that the initial time at which that trajectory leaves the detector will in general be different from the initial time for 
the unperturbed trajectory. 

In Appendix^ we derive the following useful formula for the perturbed action associated with a perturbed classical 
solution that goes from one end of an interferometer arm to the other: 

SMt')+xi{t')]+SMt')]^ r dt'L^^Hxoit')) ~ EiAtB - AtA) - r dt'H^^\xQ{t'),p^{t')), (48) 

JtA JtA 

with tB = t and tA = t — L/v, where L is the length of the arm. The trajectory xo{t') is the solution of the classical 
equation of motion associated with the unperturbed Lagrangian L(o) that goes from one end at time tA to the other 
end at time is, and poit') is the unperturbed canonical momentum associated with the trajectory Xo{t') at the time 
t' . The perturbation xi{t') vanishes at the perturbed initial and final times. The contribution from the lowest order 
term, the first term on the right-hand side of Eq. (|48|l , is identical to the time independent case and coincides with the 
result in Eq. H40() . The last term on the right-hand side of Eq. (|48() in turn coincides with the result in Eq. (|41() . Thus, 
when we add the phase differences for every stage obtained with Eq. H48|) , we recover the same result as in Sec. IIV Dl 
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except for the contributions that come from the second term on the right-hand side of Eq. H48() . Let us look at those 
more carefuUy. We are considering a given observation time at the detector location, hence Atg vanishes for the last 
stage. Moreover, the initial time perturbation At^ for any stage is always canceled by the final time perturbation 
Ats for the previous stage. Altogether we are left only with the term that corresponds to the initial time perturbation 
at the source location, Atj, for the first stage. This gives rise to a term {E/h){Ati — At[) when subtracting the total 
phase difference for the two interfering waves. On the other hand, one also needs to include the phase difference at 
the source location between the initial times ti + Ati and t[ + At[, which corresponds to the trajectories associated 
with the first and second interfering wave respectively. From the explanation provided in Sec. IIVDI it follows that 
such a phase difference is given in the present case by 

A(/. = I (ti + Ati -t[~ At'i) , (49) 

and the terms involving Ati and At[ cancel exactly those coming from the phase difference for the two interfering 
waves when computing the phase shift between them. Therefore, the result for the phase shift is entirely equivalent 
to that obtained by following the procedure based on the Hamilton- Jacobi equation employed in Secs. lIV AIIV Dl 

Chiao and Speliotopoulos obtained a different result in Refs. P, because they fixed the initial and final times 
for each stage to the same value as in the unperturbed case. Since they also considered a fixed initial velocity and 
position, the perturbed trajectory no longer coincide with the detector location at final time, but the authors did not 
realize that they were implicitly using a perturbed trajectory with the wrong boundary condition (additional details 
on this point can be found in Appendix^). As explained in the previous paragraph, when the appropriate boundary 
conditions are employed, including perturbed initial and final times for each stage, the result obtained in Sec. lIVDl is 
recovered. 



V. DISCUSSION 

We divide our discussion into two parts. First, we will examine the physical implications of the result found in 
the previous section. Next, we will discuss some fundamental aspects about the use of atom interferometers for the 
detection of gravitational waves. We will not assume any specific scheme for the interferometer {e.g. its geometry), but 
keep the discussion general by addressing the main qualitative features independent of the particular configuration. 

The phase shift between the two interfering waves in an interferometer is obtained by adding the phase difference 
(obtained in Sec. lIV|l between the ends of every arm and comparing the result for the two trajectories that correspond 
to the two interfering waves. Since the phase difference is given by an expression of the same form for all arms, the 
phase shift between the two waves will also have three contributions qualitatively similar to the first three terms in 
Eq. (|39|l . The first contribution is proportional to the difference between the total length of the paths for the two 
interfering waves and is the only term in the absence of gravitational waves. The third contribution encodes the entire 
effect of the gravitational wave on the propagation of the atoms if the mirrors were freely suspended. Finally, the 
second contribution accounts for the fact that we are considering an interferometer with rigid arms. 

One can consider two regimes depending on the ratio of the transit time of the atoms to the period of the gravitational 
wave. We will refer to the case in which the period of the gravitational wave is smaller than the transit time of the 
atoms along the arms of the interferometer as the high frequency regime. In that regime the second contribution 
mentioned above dominates. Correspondingly, the low frequency regime refers to the case in which the period of the 
gravitational wave is greater than the transit time of the atoms. In that case the second and third contributions 
become comparable, but since they carry opposite signs their effects cancel out. 

Despite the obvious differences between atom and light interferometers, e.g., the atoms in MIGO are non-relativistic, 
whereas photons in LIGO are highly relativistic, the expressions for the phase shift in the two cases are formally 
equivalent, except for the second contribution, provided that one 1) replaces the de Broglie wavelength A with the 
laser wavelength Aph and v with c (the speed of light), 2) uses the effective length for the arms of the laser interferometer 
obtained by multiplying the arm length by the number of times photons travel back and forth inside the Fabry-Perot 
cavities, b. The second contribution mentioned above is absent for LIGO because its mirrors are freely suspended. If 
its arms were rigid, the whole expression would be equivalent for atom and laser interferometers. 

Laser interferometers with freely suspended mirrors such as LIGO are mostly sensitive in the low frequency regime, 
where the expression for the phase shift becomes (5L/Aph)ft-^j^^■'?^■' times an oscillatory factor of order one which is 
a harmonic function of cot, where t is the observation time and uj is the angular frequency of the gravitational wave. 
On the other hand, for high frequency gravitational waves the phase shift becomes proportional to {c/coXphjhijn^n-' 
times shi{ujbL/c) and a harmonic function oicut. Therefore, the phase shift in the high frequency regime is reduced by 
an amount which corresponds to replacing the effective length of the arm with the distance traveled by the photons 
during one gravitational wave period, which is smaller, in that regime, than the effective length of the arms. This 
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can be qualitatively understood as follows: when several oscillations of the gravitational wave take place while each 
photon travels along the interferometer arms, the effect is averaged out and only the last cycle of the gravitational 
wave counts. 

In contrast, for most conceivable situations, MIGO would typically be implemented in the high frequency regime, 
since v <^ c for non-relativistic atoms and the transit time becomes much larger than in laser interferometers. In that 
regime, the second contribution dominates and the phase shift becomes L/X times the gravitational wave amplitude 
and a harmonic function of ujt. MIGO is not affected by that averaging-out effect for high frequencies as laser 
interferometers because its arms are rigid. Thus, it is better to consider rigid arms for MIGO rather than freely 
suspended mirrors because it works in the high frequency regime'*. On the other hand, for low frequencies the second 
and third terms become comparable but with opposite sign and cancel out. This is the same reason why LIGO would 
not work if its mirrors were not freely suspended. 

However, even though an atom interferometer with rigid arms would not exhibit the averaging-out effect mentioned 
above, the oscillatory factor limits the amount of time during which the measurement can be performed. For a fixed 
flux of atoms, that means that, as the gravitational wave frequency increases, the maximum value for the total number 
of atoms involved in the measurement decreases and shot noise becomes important, eventually limiting the sensitivity 
of the interferometer. 

Once the correct result for the phase shift, obtained by applying Eq. (|39|l to the particular configuration of interest, 
is employed, the major claims by Chiao and Speliotopoulos about the capabilities of MIGO no longer hold. First, one 
can immediately see from Eq. (39) that it is better to use fast atoms than slow ones. Second, although the averaging- 
out effect is not so drastic as for laser interferometers, the phase shift does not increase with the gravitational wave 
frequency. In fact, if one considers gravitational waves which correspond to the low frequency regime for LIGO, the 
expression for the phase shift is given by &Lligo / '^ph times the amplitude of the gravitational wave and a harmonic 
function of cot. On the other hand, that would correspond to the high frequency regime for MIGO in most conceivable 
situations, and the phase shift would be Lmigo /-^ times the amplitude of the gravitational wave and a similar harmonic 
function of ujt. Thus, under equal conditions (i.e., if one had A = Aph, ^migo = ^ligo and assumed in addition that 
the atoms could bounce back and forth the same number of times) the sensitivities of MIGO and LIGO would be 
comparable, so that the great enhancement purported in Ref. does not really exist. 

The assumptions that Xmigo — -^ligo and the atoms can bounce a similar number of times as the photons in 
the Fabry-Perot cavities are of course unrealistic. Nevertheless, one might still wonder whether a feasible range 
of parameters exists for which the sensitivity of MIGO could be comparable to that of laser interferometers. The 
main advantage of atom interferometers is that much shorter wavelengths can be used, whereas photons with much 
shorter wavelength would be difficult to work with. In theory, this could make it possible to decrease the size of 
the interferometer to more realistic values. In practice, however, there are a number of noise sources involving the 
position of the mirrors which cannot be overcome by increasing the idealized sensitivity of the interferometer, so that 
the displacement induced by the gravitational wave needs to be above a certain minimum value, which implies a 
minimum size for the interferometer. In fact, that is the reason why the size of LIGO could not be reduced simply by 
using shorter wavelength photons, since its sensitivity at low and intermediate frequencies is limited by seismic and 
suspension thermal noise. The fact that MIGO has rigid arms rather than freely suspended mirrors does not seem 
to improve the situation drastically. Finally, as should be quite obvious from the beginning, the major limitation for 
an atom interferometer is the maximum flux of atoms that can be produced as compared to the huge flux of photons 
available in a laser interferometer. Therefore, for realistic situations MIGO would likely be limited by shot noise. 

We close this section by pointing out that our conclusions are essentially the same both for interferometers which 
have a spatial projection of their arms with non- vanishing area and for interferometers with a Michelson-type configu- 
ration such as LIGO. According to Ref. [I4j . this distinction is crucial for static spacetimes with nontrivial curvature, 
because only the former type can be employed to measure the spacetime curvature. Nevertheless, it is not so im- 
portant for linearizedgravitational waves propagating on a Minkowski background. This point bears strongly on the 
claims made in Ref. [3 which depend on selecting the right expression for the phase shift among the different results 
obtained there. The authors chose the expression that coincides with the result for the static case when the static 
limit of the gravitational wave is considered. However, the static limit of linearized gravitational waves (propagating 
on a Minkowski background) is flat spacetime itself and, hence, its curvature vanishes. The static limit is trivial and 
cannot be used for that purpose, whereas for non-vanishing gravitational wave frequencies, even if they are small, the 



* In a spacetime corresponding to a harmonic linearized gravitational wave with wavelength Aqw it is only possible to have rigid objects 
with a size along the transverse direction of the order of Aqw at most. For an atom interferometer it is possible to fulfill that condition 
and still be in the high frequency regime because v <^ c. On the other hand, for a laser interferometer that condition can only be 
satisfied while being in the high frequency regime provided that the photons bounce back and forth many times inside the Fabry-Perot 
cavities. 



13 



result of Ref. |lj| for the static case is no longer valid. Anyway, when the correct procedure is employed to compute 
the phase shift, the result is the same for different coordinate systems, as one could anticipate from the principle of 
general covariance. To begin with, there is no need to find a criterion for selecting the right result. 
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APPENDIX A: ACTION FOR PERTURBED TRAJECTORIES AROUND A CLASSICAL SOLUTION 

In Sec. lIVEl we explained that one can obtain the phase difference between the two ends of an interferometer arm 
at times Ia + At^ and ts + by evaluating the action for the perturbed trajectory that satisfies the perturbed 
equation of motion and goes from one end at time tA + At a to the other end at time is + Ais. Up to first order in 
the perturbation, the action is given by 

So[xo{t') + Xi{t')] + Si[xo{t')] = dt'L(°\xo{t')+Xi{t'))+ dt'L^'\xo(t')), (Al) 

JtA + AtA JtA 

where xo{t') is the solution of the unperturbed equation of motion that goes from one end at time tA to the other end 
at time ts, and xi{t') is the perturbation of the trajectory. We have assumed that the perturbations of the initial 
and final times are of the same order as xi{t') and the Lagrangian perturbation L^^\ In this Appendix we will derive 
a useful expression for that action. In doing so, we will also point out the mistake made in Refs. '7,^. 
Let us start with the action for a general Lagrangian L 

S[x{t)] = / L{x{t))dt, (A2) 



where x{t) is a solution of the classical equation of motion associated with the Lagrangian L, and add a small 
perturbation Sx(t) to x{t). The usual variational treatment yields the following result up to linear order in Sx(t): 

SS ^pbSxb -paSxa, (A3) 

where the subindices A and B denote quantities evaluated at the initial and final time respectively. On the other 
hand, if we perturb the initial and final times, the perturbed action 

tB+AtB 

L{x{t))dt (A4) 

tA+AtA 

minus the original action, up to linear order in At a and Ai^, becomes L{x{tB))AtB ~ L{x{t a)) At a- When both kinds 
of perturbations are considered simultaneously, the final position changes, up to linear order, by Axb — Sxb +isAiB 
(considering that 5xb and isAis are of the same order). A similar result also holds for the initial time. Altogether 
the change in the action is given by: 

6S = {pbAxb - VBXBAtB + L{x{tB))AtB) - {vaAxa - PAXAAtA + L{x{tA))AtA) , (A5) 

which can be equivalently rewritten as 

SS = [pbAxb - H{xB,PB)AtB) - {paAxa - H{xA,PA)AtA) . (A6) 

Applying Eq. (|A6|) to the case in which the Lagrangian is L^*'^ and the perturbation is 5x{t) ~ xi{t), we get 

dSo = (pbAxb -paAxa) - (H^°\xB,PB)AtB - H^°\xA,PA)AtA) , (A7) 

where pA and pb are the momenta at the initial and final time for the unperturbed trajectory. Since the emission 
and observation points are fixed (when using rigid coordinates), we should take Ax a = Axb — 0. Hence, the first 
term on the right-hand side of Eq. |IA1|I can be rewritten as 

dt'L^^'^xoit')) - (H^"\xB,PB)AtB - H^"\xA,PA)AtA) - /*" dt'L(o)(a;o(i')) " E{AtB - Ai^), (A8) 

JtA ^ ' JtA 
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where we have taken into account that xq {t) is a constant-energy sohition of the unperturbed equation of motion with 
energy E. 

In addition, one can show that L^^'^{xq{x,p),xo{x,p)) = —H'^^\x,p). This can be seen as follows. Consider the 
relation between momentum and velocity for the perturbed Lagrangian^: x{x,p) — xq{x,p) + xi{x,p), where xq{x,p) 
gives the velocity in terms of the momentum when only the unperturbed Lagrangian is considered. Taking this into 
account, we have 

LW(x,i(x,p)) = L(")(x,io(x,p)) + = i(")(x,io(x,p)) (A9) 

We can then obtain the Hamiltonian: 

H{x,p) = p{xo{x,p) + xi{x,p)) - L'^^\x,x{x,p)) ~ L'^^\x,x{x,p)) 

= p{xo + x^)~L^''\x,xo{x,p))-px^-L'^^^ =H^''\x,p)-L^^\x,p), (AlO) 

and verify that H'^^\x,p) = —L'-^\xo{x,p),xo{x,p)). Together with Eq. I|A8|I this implies that Eq. (|Aip can be 
re-expressed as 

So[xo{t') + xi{t')] + S,[xo{t')] = f'\t'L^°HMt')) - EiAtB - AtA) - f\t'H^'\xoit'),po{t')). (All) 



In Refs. ll|,|9jj the authors always considered At a = AIb = 0. However, since the atoms are emitted at constant 
velocity (in rigid coordinates), in general one will have Axb ^ for the first stage, and similarly for the remaining 
stages. This violates the boundary conditions that should be satisfied when computing the phase difference between 
the ends of an interferometer arm: Axa — Axb = 0. Furthermore, when evaluating the action for the perturbed 
trajectory under this conditions, the term —E{AtB — AIa) in Eq. |IA11|I is actually replaced by pbAxb — PaAxa- 
One can explicitly check that this indeed corresponds to the extra term obtained in Refs. P, as compared to our 
result. 



APPENDIX B: POSSIBLE EFFECT DUE TO THE REFLECTION OFF THE MIRRORS 

In this Appendix we explain why, contrary to the claim in Ref . ^3 > the possible effects due to the reflection off the 
mirrors cancel out in the final result for the phase shift. 

It should be stressed that, when computing the phase difference for each arm, any effects due to the change of 
the classical trajectories as a result of the reflection off the mirror are already included in the treatment provided 
in Appendix ^ Indeed, the results there are valid for an arbitrary small perturbation of the classical trajectory 
employed to compute the phase difference. Therefore, it is clear from Eq. (|Alip that the only contribution that 
follows from integrating the unperturbed Lagrangian along the perturbed trajectory is —E{AtB — AIa)- Note that 
we can take Axa = Axb = not only for rigid coordinates but also for geodesic coordinates. The position of the 
mirror in geodesic coordinates for a rigid arm is given by Eq. (|36f) and, although the time shifts AIa and AIb can 
give rise to non- vanishing Axx and Axb, those contributions arc of higher order in hij, because At a and AtB are 
already of order /ly , and can be neglected. 

In fact, one can easily check that Eq. IjAlljl is equivalent to Eqs. (4) or (5) in Ref. However, the authors of 
Ref. ^3 did not realize that the contributions —E{AtB — At a) for each piece of the trajectory cancel out in the 
final result for the phase shift. This follows from the discussion in the next to last paragraph of Sec. IIV El For a 
fixed observation time, adding those contributions for the several pieces of the trajectory associated with one of the 
interfering waves gives EAti, where Aii is the shift of the time at which the perturbed classical trajectory leaves 
the source as compared to the unperturbed one. An analogous result, EAt[, is obtained for the other interfering 
wave. These two contributions arc, nevertheless, canceled by the contribution from the phase difference at the source 
location, which is given by Eq. (|49|l . This is schematically represented in Fig. ^ where only one of the perturbed 
trajectories is shown. In order to compute the phase shift at the observation time to, one needs to evaluate the action 
for the classical trajectory that leaves the source at time ti, is reflected off the mirror and reaches the detector at to, 
and then compare it with the result of adding the action for the trajectory that leaves at t2 to the phase difference 



^ From now on in this Appendix, all the equalities should be regarded as valid only up to linear order in the perturbation. 
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to 




FIG. 1: Schematic spacetime diagram representing the two trajectories employed to compute the phase shift in an atomic 
interferometer with a Michelson-type configuration. The vertical axis corresponds to time and the horizontal axis represents 
the positions along the two arms. Each side of the time axis represents a different direction in space which corresponds to each 
one of the arms. The continuous lines correspond to the unperturbed trajectories, and the dashed line to a perturbed one. In 
this particular example the perturbed trajectory is identical to the unperturbed one before bouncing off the mirror, but its 
velocity is higher after the reflection. 



between ti and ^2 at the source location. When a perturbed trajectory is considered (dashed line), the change in the 
action for that trajectory is exactly canceled by the phase difference between t'^ and ti at the source location. Thus, 
it is actually irrelevant whether the perturbed trajectory or the unperturbed one is used when computing the phase 
shift. 

We conclude this Appendix with several remarks. First, our previous discussion only makes sense if the difference 
between the initial times for the two trajectories is smaller than the coherence time of the wave-packets emitted by 
the source, but this is a necessary condition for the interferometer to work. Second, if the contribution due to the 
phase difference at the source location is not taken into account, one does not get the right result for the phase shift 
in an interferometer with different arm lengths even in the absence of gravitational waves (the contributions from the 
last term in Eqs. H39|) or H42|l would not be canceled by the source contribution given by Eq. (|43|l 'l. In fact, for a laser 
interferometer the phase shift obtained by considering light rays in the geometrical optics approximation would be 
entirely given by the difference of times at which the light rays leave the source. Finally, instead of considering the 
phase difference between different times at the same location for the source, one can alternatively consider classical 
trajectories with the same initial time but slightly different initial positions and include the phase difference between 
those two positions for the initial state (the wavefunction at that initial time) . This approach was employed in Ref. [l5l | 
and gives an equivalent result^. 
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